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Development of microsatellites in Labisia pumila 
(Myrsinaceae), an economically important 
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Lee Hong Tnah^'^, Chai Ting Lee^, Soon Leong Lee^, Chin Hong Ng^, 
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• Premise of the study: The exploitation of Labisia pumila for commercial demand is gradually increasing. It is therefore impor- 
tant that conservation is prioritized to ensure sustainable utilization. We developed microsatellites for L. pumila var. alata and 
evaluated their polymorphism across var. alata, var. pumila, and var. lanceolata. 

• Methods and Results: Ten polymorphic microsatellites of L. pumila were developed using the magnetic bead hybridization 
selection approach. A total of 84, 48, and 66 alleles were observed in L. pumila var. alata, var. pumila, and var. lanceolata, 
respectively. The species is likely a tetraploid, with the majority of the loci exhibiting up to four alleles per individual. 

• Conclusions: This is the first report on the development of microsatellites in L. pumila. The microsatellites will provide a good 
basis for investigating the population genetics of the species and will serve as a useful tool for DNA profiling. 
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Labisia pumila (Blume) Fern.-Vill. (Myrsinaceae) is a small 
understory shrub that is widely distributed in the tropical forests 
of Malaysia, Indonesia, Thailand, the Philippines, and Myan- 
mar (Sunarno, 2005). Eight varieties are recognized (Sunarno, 
2005), of which L. pumila var. pumila, L. pumila var. alata 
(Scheff.) Mez, and L. pumila var. lanceolata (Scheff.) Mez are 
commonly found in Malaysia. These varieties are morphologi- 
cally distinct from one another in terms of their petiole and leaf 
characteristics. Among the Malay communities, these varieties 
are collectively known as kacip fatimah, which has long been 
used as the traditional medicine for the treatment of pre- and post- 
partum complications, menstrual disorders, dysentery, rheuma- 
tism, flatulence, and gonnorrhea (Burkill, 1966; Jaganath, 2000). 

To date, the exploitation of L. pumila for commercial de- 
mand, particularly in the pharmacological and cosmeceutical 
applications, is gradually increasing. Attention to conservation 
should therefore be prioritized to ensure sustainable utilization. 
Despite the importance of L. pumila, the availability of genetic 
information for the species is still very limited. Only two ge- 
netic variability studies of the species have been reported using 
dominant markers (Bhore et al., 2009; Farah Fazwa et al., 
2013). In this study, we report the development of 10 micro- 
satellite loci in L. pumila var. alata, and we evaluate their 
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polymorphism across var. alata, var. pumila, and var. lanceo- 
lata. Microsatellites are preferred markers because of the nature 
of their codominant inheritance, high abundance, extent of al- 
lelic diversity, and the ease of assessing the size variation by 
PGR with pairs of flanking primers (Weising et al., 2005). 

METHODS AND RESULTS 

Leaf samples of 25, 20, and four individuals of L. pumila var. alata, var. 
pumila, and var. lanceolata, respectively, were obtained from Pasoh Forest Re- 
serve (2°58'N, 102°18'E). An additional six individuals of var. lanceolata were 
obtained from the Ethnobotanical Garden of the Forest Research Institute Ma- 
laysia (FRIM). The voucher specimens of these three varieties were deposited 
in FRIM Herbarium (KEP; barcode numbers 223663-223665). Total genomic 
DNA was extracted from fresh leaves of L. pumila using a modified cetyltrimeth- 
ylammonium bromide (CTAB) protocol (Murray and Thompson, 1980) and 
further purified using the High Pure PGR Template Preparation Kit (Roche 
Diagnostics GmbH, Penzberg, Germany). 

A genomic library enriched for dinucleotide CT and GT was constructed 
following the approach of Kijas et al. (1994). Approximately 5 jig of genomic 
DNA was obtained from an individual of L. pumila var. alata from the Ethno- 
botanical Garden of FRIM. After digestion with Ndell (Promega Corporation, 
Madison, Wisconsin, USA), the digested fragments were electrophoresed on 
2% agarose gels with a 100-bp DNA ladder (New England Biolabs, Ipswich, 
Massachusetts, USA). Fragments of 300-1000 bp were excised and ligated into 
Sau?>A\ cassettes (TaKaRa Bio, Otsu, Shiga, Japan). After ligation, the nicks 
were repaired using DNA polymerase I (TaKaRa Bio). The cassette-ligated 
DNA was enriched for microsatellite repeats via hybridization to 5'-biotinylated 
(CT)i5 and (GT)i5 probes and retrieved using magnetic beads coated with 
streptavidin (Promega Corporation). The selectively recovered hybrids were 
reamplified using CI cassette primers, digested with Ndell, cloned into pUCl 18 
Bam Hl/BAP vector (TaKaRa Bio), and transformed into QIAGEN EZ Com- 
petent Cells (QIAGEN Gmbh, Hilden, Germany). Insert-containing clones 
were selected by blue/white screening on Luria-Bertani (LB) agar plates con- 
taining 100 jig/mL ampicillin, 50 jiM isopropyl-p-D-l-thiogalactopyranoside 
(IPTG), and 80 |ig/mL X-gal. Plasmid DNAs of a total of 608 clones were 
amplified using the Illustra TempliPhi Amplification Kit (GE Healthcare, 
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Table 1. Description of 10 polymorphic and three monomorphic microsatellites screened in Labisia pumila van alata. 



Locus 


Repeat motif 




Primer sequence (5'-30 


Lpu02 


(GT)2o 


F: 
R: 


GCGGAAAGAAGGTATAGTGTG 
AAATATAATAGGCGACACGAGTTA 


LpuOSa 


(TC)i5 


F: 
R: 


ACAAATCCAAATTCGCACCT 
AAT GAAAAAT AAT AAT AAAC C C AA 


LpuOSb 


(TC)7 


F: 
R: 


CTTCCTGTCCTTGTTGGTTAGAGT 
AGAGGGAAAGTCGAAGTTAAT 


Lpul3 


(TG)ii 


F: 
R: 


CGAACTGAGAATATCCCCTTGCTA 
TAGTTGATGTGAGTTGGCGTAAA 


Lpul5 


(TTTTC)3 


F: 
R: 


GGGCACTTGTCTTGTAGGTTAGAG 
TCAGCGACCGAAACATCACTACTC 


Lpul6a 


(CAC)5 


F: 
R: 


GCCGAGCTACCAAGGTTAC 
GCCTCTGAGGAAAGACAAGCATAA 


Lpul6b 


(GTAT)io 


F: 
R: 


TTTCCTCAGAGGCTGGTATGTAT 
TTTCAAGGGTGTCAATCAATTAAC 


Lpu21a 


(AG)v 


F: 
R: 


AGAGGGAAAGTCGAAGTTAAT 
CTTCCTGTCCTTGTTGGTTAGAGT 


Lpu21b 


(GA)i5 


F: 
R: 


AGGAAGAAGCTGCCGTGAGA 
GGCGAAATCACCACCACAGTTAGA 


Lpu38 


(CA)8 


F: 
R: 


TTCCTACTCACTGCACATAGTTGC 
GCTGTAGCGTGTGCGAGTG 


Lpu06* 


(CA)6 


F: 
R: 


TAAAGCACCACAATTCACGCACAC 
CGCTTAGGATTTGAAGTGGGTAA 


Lpu20* 


(CT)6 


F: 
R: 


CATGCCCATACACCGGAAAT 
GAGGCGAAAAGTCACATAAACGTA 


Lpu24* 


(GGAATT)3 


F: 
R: 


CTATTATTGTGTTGATGGATTGAG 
GC C AAT T T C AT AAC C AAT C AT A 



Note: = annealing temperature. 
* Monomorphic microsatellites. 



Piscataway, New Jersey, USA) and sequenced using BigDye Terminator Se- 
quencing Kit version 3.1 (Applied Biosystems, Foster City, California, USA) 
on an ABI 3130x1 Genetic Analyzer (Applied Biosystems). 

Sequence redundancy was checked using CD-HIT (Li and Godzik, 2006), 
and microsatellite-containing sequences were then identified from the nonre- 
dundant sequences using the MIcroSAtellite Identification Tool (MISA) (Thiel 
et al., 2003). A total of 56 sequences harbored at least one microsatelHte motif. 
Forty primer pairs of 19-24 nucleotides in length were designed using OLIGO 
6.67 (Molecular Biology Insights, Cascade, Colorado, USA), based on the cri- 
teria of GC content (40-60%), melting temperature (52-58°C), and the lack of 
secondary structures. Four unrelated individuals of L. pumila var. alata were 
used for initial screening of these primer pairs. The PCR amplifications were 
carried out in 10-jiL volumes containing 10 ng of template DNA, Ix GoTaq 
Flexi Buffer, 1.5 mM MgClz, 0.3 |iM of each primer, 0.2 mM of each dNTP, 
and 0.5 unit of GoTaq Flexi DNA polymerase (Promega Corporation). The re- 
action mixture was amplified using a GeneAmp PCR System 9700 (Applied 
Biosystems) with the following cycling conditions: 4 min at 94°C; 35 cycles of 
94°C for 1 min, 50°C for 30 s, and 72°C for 30 s; followed by a final extension 
of 30 min at 72°C. Twenty-nine primer pairs that showed specific amplification 
products of expected fragment size were selected for fluorescent labeling at the 
5'-end of the forward primers with either 6-FAM or HEX. These primers were 
further screened using 25 samples of L. pumila var. alata using the same cy- 
cling conditions as above. Fragment analysis was conducted using an ABI 
3130x1 with GeneScan 400 ROX (Applied Biosystems) as the internal size 
standard. Allele sizes were assigned using GeneMapper version 4.0 (AppHed 
Biosystems). Characteristics of the microsatellites developed for L. pumila var. 
alata are shown in Table 1 . Thirteen loci yielded consistent and scorable geno- 
types; of these, 10 were polymorphic and the remaining three were monomor- 
phic. Because the majority of the loci evaluated (Lpu02, LpuOSa, LpuOSb, 
Lpul3, Lpul6b, Lpu21a, and Lpu21b) exhibited up to four alleles per individ- 
ual (Fig. 1), our results suggest that the species is likely a tetraploid. 

To assess the degree of polymorphism across L. pumila var. alata, var. 
pumila, and var. lanceolata, the 10 polymorphic microsatellites were further 
checked for variability in 55 individuals of the three varieties using the same 
PCR and genotyping protocols as above. Observed heterozygosity iH^) of each 
variety was calculated based on the allele phenotypes, following Bever and 
Felber (1992), with values being one minus the probabihty that any two alleles 
drawn at random were identical by descent (AAAA = 0, AAAB = 0.50, AABB = 
0.667, AABC = 0.S33, and ABCD = 1). Minimum and maximum bounds of //„ 



7;(°c) 


Allele size range (bp) 


Fluorescent label 


GenBank accessi 


50 


197-235 


HEX 


KF31S311 


50 


103-132 


HEX 


KF31S312 


50 


17S-211 


HEX 


KF31S312 


50 


303-342 


HEX 


KF31S313 


50 


19S-213 


6-FAM 


KF31S314 


50 


93-107 


6-FAM 


KF31S315 


50 


69-112 


6-FAM 


KF31S315 


50 


175-207 


6-FAM 


KF31S316 


50 


3SS-415 


6-FAM 


KF31S316 


50 


75-S3 


HEX 


KF31S317 


50 


14S 


HEX 


KF31S31S 


50 


315 


6-FAM 


KF31S319 


50 


159 


6-FAM 


KF31S320 



were calculated for bialleUc phenotypes (AAAB and AABB) (De Walt et al., 
2011). Expected heterozygosity (H^ was calculated using ATETRA 1.2 (Van 
Puyvelde et al., 2010), with 10,000 Monte Carlo simulations. A total of S4, 4S, 
and 66 alleles were observed in L. pumila var. alata, var. pumila, and var. lan- 
ceolata, respectively (Table 2). ranged from 0.039 to 0.S57, 0.000 to 0.793, 
and 0.000 to 0.S74 in L. pumila var. alata, var. pumila, and var. lanceolata, 
respectively. Notably, locus Lpul5 was found to be monomorphic in varieties 
pumila and lanceolata. 



I Lpu02 

2Q7 m m 225 231^ 

3000] ^ ' * ^ ' ' * ■ ^ 




al 179 al 189 al 197 a 1 203 

Fig. 1 . Electropherogram showing four alleles amplified from an indi- 
vidual at loci Lpu02 and Lpu21a. 
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Table 2. Genetic properties of 10 microsatellites of Labisia pumila across varieties alata, pumila, and lanceolata. 



Locus 




L. pumila var. alata (n = 


25) 




L. pumila var. pumila (n 


= 20) 




L. pumila var. lanceolata (n = 


10) 


A 


Hq (min) 


1 Q V^lllClW J 




A 


Hq (min) 






A 


(min) 


H (max) 




i^puuz 


1 1 


0.940 


0.953 


U.OJ / 


J 


0.783 


0.808 


U. / jo 


Q 


0.483 


0.584 


U.oUO 


LpuOSa 


9 


0.880 


0.880 


0.793 


7 


0.833 


0.833 


0.793 


8 


0.717 


0.783 


0.815 


LpuOSb 


11 


0.940 


0.940 


0.834 


6 


0.833 


0.833 


0.763 


8 


0.783 


0.850 


0.794 


Lpul3 


10 


0.706 


0.753 


0.848 


6 


0.533 


0.684 


0.695 


7 


0.370 


0.426 


0.652 


Lpul5 


2 


0.020 


0.027 


0.039 


1 


0.000 


0.000 


0.000 


1 


0.000 


0.000 


0.000 


Lpul6a 


5 


0.280 


0.340 


0.412 


2 


0.075 


0.100 


0.138 


2 


0.150 


0.200 


0.253 


Lpul6b 


10 


0.973 


0.973 


0.853 


5 


0.833 


0.833 


0.729 


9 


0.816 


0.833 


0.852 


Lpu21a 


11 


0.940 


0.940 


0.834 


6 


0.833 


0.833 


0.763 


8 


0.733 


0.817 


0.785 


Lpu21b 


10 


0.972 


0.972 


0.843 


7 


0.833 


0.833 


0.777 


11 


0.917 


0.917 


0.874 


Lpu38 


5 


0.806 


0.820 


0.731 


3 


0.125 


0.167 


0.405 


3 


0.250 


0.334 


0.591 



Note: A = number of alleles; = expected heterozygosity; (min) = minimum observed heterozygosity; (max) = maximum observed heterozygosity; 
n = number of individuals. 



CONCLUSIONS 

This is the first report on the development of microsatelhtes 
in L. pumila. The observed levels of polymorphism and genetic 
diversity suggest that, apart from monomorphic loci, these mi- 
crosatellites can serve as useful tools for DNA profiling and 
population genetic studies of the species. 
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